Quasi-one-dimensional scattering in a discrete model 
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We study quasi-one-dimensional scattering of one and two particles with short-range interactions 
on a discrete lattice model in two dimensions. One of the directions is tightly confined by an arbitrary 
trapping potential. We obtain the coUisional properties of these systems both at finite and zero Bloch 
quasi-momenta, considering as well finite sizes and transversal traps that support a continuum of 
states. This is made straightforward by using the exact ansatz for the quasi-one-dimensional states 
from the beginning. In the more interesting case of genuine two-particle scattering, we find that 
more than one confinement-induced resonance appear due to the non-separability of the center-of- 
mass and relative coordinates on the lattice. This is done by solving its corresponding Lippmann- 
Schwinger-like equation. We characterize the effective one-dimensional interaction and compare it 
with a model that includes only the effect of the dominant, broadest resonance, which amounts to 
a single-pole approximation for the interaction coupling constant. 

PACS numbers: 03.65.Nk, 34.50.-s, 32.80.Pj, 



I. INTRODUCTION 

Low-dimensional systems have been among the dreams 
of theorists for decades |1|, since they provide a per- 
fect playground where many analytical techniques be- 
come exact for a variety of models. For instance, we have 
the Bethe ansatz for the one-dimensional Bose gas with 
contact interactions [3], or bosonization for Luttinger's 
model [3|, Ij] . Simple product ground-state wave functions 
are exact in supersymmetric one-dimensional quantum 
systems, with Sutherland's model @ their main repre- 
sentative. In addition, numerical methods such as the 
density-matrix renormalization group [6] and its recent 
formulation in terms of matrix product states ill can 
handle non-integrable models efficiently and accurately 
in one dimension. 

With the advent of ultracold atom and molecule 
physics [8|, it is now possible to engineer effective two- 
and one-dimensional many-body systems by strong con- 
finement in one or two dimensions. This possibility, com- 
bined with the great degree of control of two-body inter- 
actions thanks to magnetic Feshbach resonances pj|-[ll|. 
has led to major experimental achievements. For ex- 
ample, a Tonks- Girardeau gas |I2| - a one-dimensional 
gas of impenetrable bosons - has been experimentally 
realized in an optical lattice [IJ]. In 1998, Olshanii de- 
rived how short-range interacting bosons effectively in- 
teract in the quasi-one-dimensional regime when tight 
confinement is applied in two of the space directions [14| , 
and predicted the existence of a so-called confinement- 
induced resonance under experimentally relevant condi- 
tions. This type of resonances was used in [15] to create 
and characterize a super Tonks- Girardeau gas, that is, a 
highly excited Bose gas state where the interactions are 
strongly attractive. 

Confinement-induced resonances have been recently 
studied experimentaly in detail in 16j]. There, Haller 
and co-workers found that, as the transversal trap was 



made anisotropic, a splitting and shift of an inelastic res- 
onance - only atom losses were measured - occured, and 
that additional resonances appeared as anisotropy was 
further increased. There have been various theoretical 
attempts to explain the splitting of the resonance. In 
[l7|, Peng et al. studied two-body s-wave collisions un- 
der anisotropic harmonic confinement, and found that 
in this model a shift in Olshanii's resonance does occur, 
but no splitting was observed. Moreover, the resonance 
positions are shifted from the experimental data. This 
issue has been beautifully resolved by Sala and collab- 
orators in [1^, where they showed that the splitting is 
due to non-separability of the center-of-mass and relative 
motion in two-body scattering, since in the actual exper- 
iment an optical lattice - intrinsically anharmonic - was 
used. Sala et a/.'s theory and results have subsequently 
been corroborated in [19|] . 



In this paper, we study one- and two-body scattering 
with zero-range interactions in the quasi-one-dimensional 
regime. The particles live in two spatial dimensions, 
for simplicity. To avoid unnecessary complications with 
the singular nature of the zero-range potential, we fo- 
cus on a lattice model. An interesting by-product of 
our discrete model is that the center-of-mass and rel- 
ative motion in the two-body problem are not separa- 
ble. We are therefore able to confirm the appearance of 
several confinement-induced resonances and we also pre- 
dict the occurrence of confinement-induced cancellation 
of the effective interaction. We deal here as well with 
some general situations that have not been considered 
so far: transversal confinement supporting a continuum 
of states and finite quasi-one-dimensional systems sub- 
jected to periodic boundary conditions. The latter case 
is relevant to ongoing research on ultracold atoms in ring 
geometries |20 |. 



II. POTENTIAL SCATTERING 

We begin by studying a single particle colliding with a 
zero-range potential barrier or well. This is equivalent to 
the two-body problem when it is separable into center- 
of-mass and relative coordinates, the case first studied by 
Olshanii [l4| . We first derive the relevant expressions for 
the scattering length and phase-shifts in the case of an 
infinite system, and we then obtain our results for finite 
systems and transversal traps supporting a continuum. 



A. Infinite quasi-one-dimensional space in a trap 

We study the problem of a single particle in two spa- 
tial dimensions, where only one of the directions (y) is 
trapped, while the other {x) is free. At {x,y) = (0,0) 
there is a zero-range potential of strength U. We use 
a lattice model instead of the continuum one to avoid 
unnecessary complications with the irregularity of the 
delta potential in dimensions higher than one. Moreover, 
quasi-lD physics has not been studied so far on a lattice 
and shows interesting phenomena for the non-separable 
two-body case, as we will see. The Hamiltonian of the 
system reads 



H — Hx + Hy + U5x,o6yfi, 
where (5a t is a kronecker delta and 



(1) 



Hx = ~JY,^\x + 1){x\ + \x){x+1\), (2) 

X 

Hy = -jY,{\y+l){y\ + \y){y+l\ + V {y) \y) {y\) , 



with J > the single-particle tunneling rate, x and y in- 
teger numbers (we set the lattice spacing d=l) denoting 
the position of the particle, and V{y) is a trapping poten- 
tial supporting no continuum states. The trap V{y) can 
have any shape and strength, and we begin our analysis 
with a harmonic potential, V(y) = i^y^', the properties 
of the spectrum and eigenfunctions of Hy, {En}n>o and 
{V'n}n>o, have been studied in |2ll.l2a|. 

We start by finding the one-dimensional scattering 
length a of the system. That is, we obtain the lowest- 
energy scattering solution to the stationary Schrodinger 
equation iJ^ — E'^, with eigenenergy E — —2J + Eo. 
The simplest way to proceed is by proposing the following 
exact ansatz 



^(x,y) ^{\x\~a)Mv) + 



^fe„alf'V'«(y), 

n=l 



(3) 



where |a„| < 1 must be satisfied for each n, provided 
E > —2 J. At higher quasi-momenta, |a„| may attain 
unit value for low n. In Eq. ([3]), &„ are some expan- 
sion coefficients to be determined. Two remarks on the 
form of Eq. ^ are in order, (i) The first term is the 



lowest-energy scattering solution (the analog to the zero- 
energy solution in the continuum) in one dimension [23| 
multiplied by the transversal ground state; it represents 
the quasi-lD solution we are looking for. (ii) Every term 
in the infinite sum corresponds to a virtual excitation to 
the n-th trapped transversal state., (iii) The same ansatz 
works for the continuum theory of 'l4', of course. 

We introduce the ansatz ([3]) in the Schrodinger equa- 
tion H'ii = E'^ and, for a; 7^ 0, it is satisfied provided 



J- 



.1 



En — E, 



(4) 



for all n. The above relation expresses the fact that the 
energy of the virtual bound state [23| lowers the total 
energy to E with respect to £"„. This point is not spe- 
cific to our model, but appears in general multi-channel 
scattering problems involving an energy gap to coupled 
excitations [2J]. We now use Eqs. ^ and ^ in the 
Schrodinger equation at a: = 0, and project onto the 
transversal eigenstates i/)„ [ij, [T^l , obtaining 



*(0,0) 



2J 



bn^U 



uroioy 

V4 (0)^/(0,0) 

2 Jan + E-En 



(5) 
(6) 



Inserting the relations ^ into ^1^(0,0), Eq. (|3]), and 
equating it to Eq. ([5]) we find the scattering length 



2J 



u\Mo)\ 



l-U 



\MoW 



, 2Ja„ 

n—l 



E — En 



(7) 



At low energies, we can map our problem into an ef- 
fective one-dimensional model. To do so, recall that 
the scattering length for a particle colliding with a ID 
delta potential of strength [/id on a lattice is given by 
am = -2J/Um [25, 26]. Therefore, the effective ID 
potential strength C/id is given by 



Um = 



1 - c//c/( 



CIR 



where C/ciR is the value of U at which there 
confinement- induced resonance (CIR), that is, Uyd - 
It is given by 



1/t/, 



CIR 



S 2^"" 



l^n(0)| 



E — En 



(8) 

is a 

> 00. 



(9) 



In Fig. [T]we show the effective ID interaction strength 
UiY) as a function of the bare coupling strength U/ J for a 
shallow harmonic trap with Vt/ J = 10"'^. Note that the 
expression obtained for C/id, Eq. ([8]), coincides with its 
continuum counterpart [ij] . They differ only in the spe- 
cific, model-dependent value of t/ciR- The two limiting 
cases of weak and strong U are also analogous: C/id = 
for [7 = 0, while [/id -^ -C^cir|V'o(0)P when U -^ 00. 
The major difference is that the continuum CIR is found 



to occur at a positive (renormalized) coupling strength 
(in 3D), while in our case it occurs at negative values 
(see Fig. [T]); this has to be so, since the bare coupling 
constant in the continuum is brought to the opposite sign 
after renormalization. 

We now consider the case of a finite incident momen- 
tum fc (0 < |A:| < tt). The ansatz for the lowest energy 
solution is readily generalized by simply changing |a::| — a 
into cos{k\x\+Sk) in dS]), where 6k is the phase shift to be 
determined. In Eq. ((H), the energy E is substituted by 
E{k) = — 2Jcos(fc) + Eq. After analogous manipulations 
to those carried out for the scattering length above, we 
obtain 



tan Sh 



Ucscjk) \M0)\^ 
2 J 1 - U/Ucm{k) ' 



where 



l/Ucmik) = J2 



\M0)\' 



2 Ja„ + E{k) - En 



(10) 



(11) 



The above expression reduces to Eq. © for A; = 0, 
so UciR = f^CiR(O). For finite momenta, by compari- 
son with the pure ID solution [25l . l26l |. we obtain the 
momentum-dependent effective interaction 



Uiu{k) 



u\M0)\' 

l-U/Ucm{ky 



(12) 



There is a relevant remark to make about the above 
results. First, it is clear that, while at weak potential 
strength U the effective ID properties are dominated 
by the bare coupling, in the hard-core limit these are 
dictated by the virtual excitations in the trap, that is, 
the position of the confinement- induced resonance, C/ciR- 
This means that many-body physics in the limit of hard- 
core interparticle interactions will be solely determined 
by the virtual transitions, whenever the densities are low 
enough for having all |a„| < 1 in Eq. ^ at zero and 
finite relative quasi- moment a. 



B. Infinite quasi-one-dimensional space in a trap 
supporting a continuum 

So far, we have dealt with the case of a transversal 
trapping potential that supports only bound states. Such 
a situation is obviously ideal, and we here address this 
issue within the potential scattering approach. 

We consider a single-particle described by Hamiltonian 
([T]), with H^ and Hy defined by Eq. ©, and V{y) a 
transversal potential well such that Hy supports at least 
one and at most a finite number of bound states. For 
convenience, we assume that V is symmetric {V{y) — 
V{—y)) and monotonic (l^(|yl + 1) > V{\y\)), and we 
add an energy off-set so as to render it positive, V{y) = 
V{y) + Vo, with Va = -V{0) > 0. We will consider 
V to have an arbitrary but finite range -R > 0, that is, 
V{\y\ >R) = Vo. 




FIG. 1. Inverse tangent of effective interaction as a function 
of the bare single-particle potential strength for Q/J = 10^'^. 



The symmetric bound states of the transversal Hamil- 
tonian Hy are readily obtained by calculating the roots 
of a certain polynomial of degree 2R -I- 1 if i? > 1 [2^ 
or of degree 2 for the zero-range case |23l . |25| . The an- 
tisymmetric states do not contribute to the scattering 
length/amplitude. The scattering states, with their cor- 
responding phase-shifts are also easily calculated [23| . 

The effective quasi-lD scattering lengths and phase 
shifts have the same expressions as for the trapping case, 
Eqs. ([7]) and ([TOl) . while special care has to be taken 
in the sums of Eqs. ^ and (fTTj) . First of all, no- 
tice that if Hy had a purely continuous spectrum, the 
quasi- ID picture would break down and the sum in Eq. 
PT|) would not be well-defined; in such case genuine two- 
dimensional scattering solutions must be considered us- 
ing, e.g. the Lippmann-Schwinger equation in quasi- 
momentum space. Fortunately, this situation happens 
only if y = 0, and we need not worry about repulsively- 
bound states |25l - l27l | above the transversal continuum if 
we choose V{R) < Vq. 

The symmetric scattering states of Hy have the form 



Aiy) 






+ e,)ii\y\ 

\y\<R. 



> R, 



(13) 



In the above equation, q is a quasi-momentum, dq are 
the phase shifts, (pq is the inner part of the scattering 
wave function to be determined, and iV is a quantization 
length unit which ensures normalization for N ^ R. 

In order to be able to perform sums over continuum 
states, we need their density of states. Applying the 
boundary condition ipq{L) = 0, with 2L + 1 = A^, to Eq. 
P^ . we obtain the usual discrete set of allowed momenta 
Q ~ qn — (2« + l)7r/2L — 9q/L. Following a procedure 
parallel to that of ref. |28| (see chapter 10), we obtain 
the density of states g{q), 



gii) 



L 

■K 



TT dq 



(14) 



We express the discrete sum of any function F as 

J2 FiQ) = E ^(9)5(g)Ag ^ r dqFiq)g{q), (15) 

q q •''^ 

where we have used that g{q)Aq — 1 and we have as- 
sumed that the limit L —>■ oo can be taken to replace 
the sum by an integral over the Brillouin zone. We then 
obtain the relevant integral over continuum states 



sik)^Yl 



\MoW 



E{k) - E{q) + {Eo - Vo) + 2Jaq 



i- I dq- 



IM^W 



2n J_, ^ E{k) - E{q) + (^0 - K,) + 2 Ja, 

1 r . i0,(o)p50,/ag 



lini 



dq- 



L^oo 2ttL J_^ E{k) - E{q) + {Eo - Vo) + 2 Ja, 

(16) 

In the above equation, E{k) = — 2Jcos(fc) + Eq is 
the energy of the target quasi-lD state (fc is the one- 
dimensional quasi- momentum), and aq satisfy Eq. (j4]) 
with their corresponding energies changed as En — >■ 
Eiq) + Vq - Eq. The last term in Eq. ^ is kept 
since it may contribute in the case of a sharp resonance 
i,dq0q = oo). If a sharp resonance at q = qp contributes 
to dni, then dq0q oc LS{q - qo) for large L [H. 

We now have all we need for the solution to our prob- 
lem, so that the position C/cm of the CIR is given by 



M-l 



l/Ucmik) = E 



l^n(0)| 



2 Ja„ + E{k) - En 



+ Sik), (17) 



with S given by Eq. (jl6|) . and A4 the number of bound 
states supported by the trap. 

To illustrate how continuum states of Hy are included, 
we consider the simplest possible transversal potential 
^{y) = yo — Vo^yfi, with Vb > 0, for which the above ex- 
pressions can be obtained analytically. In this case, the 
number of bound states is A^ = 1. We calculate the ef- 
fective quasi-lD scattering length, so we set E = E{0) — 
-2 J + Eo, and we have [M HI ^o = -y/vf+Ij^+Vo- 
The phase shifts of transversal scattering states are given 
by ta.n6q = Vo/{2Jsinq), and therefore there are no 
sharp resonances. The values of aq are given by 



^f{q,Vo/J)-Vfi<l.Vo/J)^-l, (18) 



with f{q, r) = 1 — cos q + ^(r/2)2 + 1, and hence satisfy 
< a, < 1 for all q for Vq > 0. From Eq. ([16|) we then 
have 



I r 1 

5(0) - - / dq 

27r J_^ 1 + J-: 



{2 J sin q)'-' 



2J{cosq-l) + 2Jaq 

(19) 

which is finite, as we expected. Note that, in this very 
simple model, it holds that l/C/cm = 5(0). In Fig. [2] 
we show the resonance position C/ciR ^-s a function of the 
model trap strength Vq. 
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FIG. 2. Position of the CIR as a function of the trap depth 
for a model trap supporting a single bound state and a con- 
tinuum (see text). In the inset, an augmented portion of the 
figure shows the resonant positions before the transversal trap 
reaches the extreme strong-coupUng limit. 



C. Finite quasi-one-dimensional space in a trap 

We here construct quasi- ID scattering states with pe- 
riodic boundary conditions (PBC) in the x direction. 
This represents the case of a particle in a ring trap of 
circumference L on the surface of a very long cylinder 
wrapped around the y-axis, interacting with a zero-range 
potential at (x, y) = (0,0). The following ansate repre- 
sents the exact quasi- ID scattering states (with energy 
E = -2Jcos(fc) -I- ^o) subjected to PBC 



*(a;,2/) = cos{k{\x\ - L/2))Mv) 



(20) 



OO 

E K [aL"l + a. 



n ^n 



V'«(2/), 



where x is defined mod L, 6„ are expansion coefficients, 
and a„ satisfy the relations —J{an + l/o^n) -I- En = 
— 2Jcos(fc) -I- i?o, with |a„| < 1 for all n. Note that 
in this case no scattering length can be defined, since for 
any non-zero interaction strength U and L < oo, fc = 
is not an allowed quasi-momentum. We remark that the 
phase shift has been fixed so that the state satisfies PBC, 
and is given by the same expression as that in the Lieb- 
Liniger model [4] particularized to a lattice; note that we 
need one further condition to determine the value of fc, 
given below. 

After introducing ^ into the Schrodinger equation 
H"^ = E"^, with E = -2Jcos{k) + Eo, and following a 
procedure parallel to the previous subsections we obtain 
a nonlinear equation for the ground-state momentum k 



2Jsm{k)tan{kL/2) = 



(21) 



^ + uj:Zi 



H'AOW 



{E„~E){l+al)-2J{a„+af) 



Allowed quasi-lD energies (with the zero mode of the 
trap Eq as a reference) for different system sizes are plot- 
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FIG. 3. Allowed scattering energies — 2Jcos(A;) as functions 
of the single-particle interaction strength f// J for a finite ID 
space with periodic boundary conditions, of length L = 10 
(black solid lines), L = 50 (red thick dashed lines) and L = 
1000 (blue thin dashed lines), for fi/J = 10"^. 



ted in Fig. |3l where finite-size effects are apparent. A 
CIR occurs every time the energy crosses the "fermion- 
ized" energy, that is, the energy of a non-interacting, 
antisymmetric quasi-lD state {k = (2n + l)7r/L), which 
can only happen for [/ < 0. 

A continuum of transversal states is again dealt with 
by Eqs. (fT6|) and p?|) . where we now have to replace the 



terms 2Ja„ and 2Jaq by 2J(q;„ -I- a^ ^) and IJiotq 
a^~^), respectively. 



III. TWO-BODY SCATTERING 

In the system we are considering in this article, two- 
particle collisions cannot be studied via separation of the 
CM and relative coordinates. It is therefore of interest 
to obtain the quasi- ID properties of two-body scatter- 
ing, for which non-separability effects may be relevant. 
These have proved to be fundamental [iSl to explain the 
experimental results of Haller et al. ^^ . 

We consider the quasi-lD space to be an infinitely long 
line, with the dynamics governed by the following Hamil- 
tonian 



H — Hr 



H^ 



H,, 



Hni 



U dx-i.x2"yi,y2^ l^^j 



where {xi,yi), i = 1,2, represents the position of particle 
i in two dimensions, H^i, Hy. are as defined in Eq. ([2]), 
and U is the two-body interaction strength. 

We solve the stationary Schrodinger equation H'^ = 
i?^, with ^ corresponding to a quasi-lD scattering state. 
The exact symmetric ansatz for two bosons (or two spin- 
singlet fermions) now has the form 

^{xi,yi;x2,y2) = e'-'^-^ cos(fc|a;i - X2\ + 4)'0o(yi)^o(y2) 



where X = {xi +X2)/'2., K is the (conserved) total quasi- 
momentum in the a:-direction, and where J- represents 
the closed multi-channel scattering part of the wave func- 
tion, 

J^{xi,yi;X2,y2)^ ^ bnun2a-ni.n2''^S[lljni1pn2]{yi,y2)- 
ni<n2 

(24) 
In the above, 5 is a symmetrization operator [29| . and 
the sum excludes (ni, 712) = (0, 0). After introducing (^5)) 
into the Schrodinger equation, we obtain the condition 
for a„j_„2 so that the ansatz is asymptotically correct, 

1 + "n 
— Jk '-^ — E — Eni — En2, (25) 

with Jk = —2Jcos{K/2) the collective tunneling rate 
[25I |26| . Following a methodology analogous to that for 
potential scattering, we obtain the expansion coefficients 



Uln,., 



E + 2JKani,n2 ^ £-«! — -E'n2 



where the (unknown) functions I are defined as 

00 
2'ni,n2 = X! '5'[V'*iC2](y,y)*(0,0;y,2/). 

y=-oo 

The phase shift is given by the equation 
2Jk sinksinSk — — C/Xq.o- 



(26) 



(27) 



(28) 



Defining now the overlaps between different non- 
interacting states in the trapped direction as 

00 
7^(r^l,n2;TOl,r7^2) = ^ S[ip:;^^ip^^]{y,y)S[ipm,'il;m2]{y,y), 

y— — oo 

(29) 
the final non-linear system of equations in the unknowns 
1ni,n2 is found to be 



X.., 



7e(ni,n2;0,0)Wl 



/ C/Io., 






\2Jk sin/c 
7^(ni,n2;mi,m2) 



H/ -f- ^J Kdtmi^mQ, -^nii ^1712 






(30) 



where the sum is restricted, as before, to (7711,7712) ^ 
(0,0). 

At very low momenta, the effective ID scattering 
length dominates the physics and, to calculate it, we only 
have to replace cos(fc[x| + 6) hy \x\ — a in Eq. ([23| . The 
system of equations we obtain now reads 



e' T{xi,yi]X2,y2), 



(23) 



ifii,n2 = 7?.(77i,7l2;0,0) 



^ E 



(31) 



7?.(7il, 772; 7771, 7772) 



m\<m2 



^ y ^^ K^mi,m2 ^mi ^1712 



All quantities in the above equation have the same mean- 
ing as for the non-zero quasi-momentum case, with the 
sum excluding (toi,TO2) — (0,0), and the relevant com- 
ponent being now Iq.o = Uiy)/U — 2JK/aU. Note 
that Eq. (|3T|) is a system of linear equations, in con- 
trast with Eq. (pOl) . Eq. ([31]) resembles the Lippmann- 
Schwinger (LS) equation- although with no singularities 
in the "integrand" - if the following identifications are 
made: (i) (ni,n2) = Qn serves as a discrete momentum 
state, (ii) C/7?.(q„;qm) = V{(\n](\m) is an effective po- 
tential, (iii) 1/{E + 2Jxaq^, — E^^^ ~ ^ma) pl^ys the 
role of a non-interacting Green's function G", and (iv) 
Ulq^^ = T(q„,0) represents the T- matrix. We can for- 
mally write Eq. dH]) as T = V + VG°T; in particular, 
this analogy allows us to write its Born series as 



T = V + VG°V + VG^VG^V 



(32) 



which is ensured to converge whenever U is smaller in 
magnitude than the first (smallest in magnitude) CIR 
coupling t/^jj^. 

We will see that, due to the non-separability of the 
CM and relative coordinates in the trap channels, there 
are more than one CIR and, interestingly, there are true 
quasi-lD resonances (a — ;> oo). In order to understand 
these, it is instructive to first study an easier problem in 
which the trap has only two states. This can be mod- 
elled by choosing V{y) — 2V5y^i - a discrete, transversal 
step function ~ with open boundary conditions in the y- 
direction, '0n(y = —1) = V'«(2/ = 2) = 0. With this 
choice, both the effective Green's function and poten- 
tial in the LS equation T = V + VG^T can be calcu- 
lated analytically, and its numerical solution is trivial. 
In Fig. m we show Uid as a function of U and the po- 
sition of the confinement-induced resonances for a fixed 
value of V/J = 1. There, we see that two CIRs and 
one true quasi-lD resonance are present. The rightmost 
resonance &t U ~ C^ciR ^^ very sharp, while the other 
resonance aXU = Uq^yi is wide and its shape is similar to 
that of potential scattering or, for what matters, of the 
case of separable CM and relative coordinates. Based on 
the knowledge gained in Sect. [TTl we may expect that, 
approximately, the position of the leftmost resonance is 



given by U^ 



(2) 
CIR 



lim^^oo t^i_D/7?.(0, 0,0,0), since the 



jump in [/id from the right to the left (or viceversa) of 
^ciR i^ rather small. To test this hypothesis, we start 
from the analogous expression to Eq. ([S]) for the two- 
body case, which we may call single-pole approximation 
(SPA), 



U- 



ID 



[/7e(o, 0,0,0) 

1 - u/u, 



(2) 
CIR 



(33) 



We expand the above expression for U ^ t^ciR? a-nd 
perform a least squares fit Uyd k, ci + C2/U to the nu- 
merically calculated effective interaction. Equating the 




FIG. 4. Inverse tangent of effective interaction as function 
of the bare two-body interacting strength for the two-state 
trap model with V = J . Red dots represent exact results, 
blue solid line corresponds to single-pole approximation (see 
text), and dashed vertical lines mark the CIR positions. Inset: 
zoomed- in region near f/ciR' with axes labels the same as for 
main figure. 



SPA (l33l) to zero-th and first orders to the fit, we obtain 



U. 



(2) 
CIR 



Cl 



i?(0, 0,0,0) 



^C2 



i?(0, 0,0,0)' 



(34) 



Since the SPA of Eq. (|33|) is not exact, the last equal- 
ity in the above equation is not satisfied, but the differ- 
ence between both sides of the equality provides a good 
test for the validity of the SPA, and a qualitative er- 
ror estimation. We have done a least squares fit from 
U/J = -1000 to U/J = -900, using 50 equidistant 
points; the trapping strength is set to V/J = 1. The 
results we obtain are Ci = 21.57J, C2 — — 661.22J^, and 
we have i?(0, 0,0,0) = 3/4. With ci (cs), from Eq. ^ 

we obtain a value of U^cm/ J ~ -28.76 (-29.69). Qual- 
itatively, we locate the resonance position at « — 29.2J, 
with an error of 0(J). It compares favorably with the 

(2) 



exact location [/qjj^ = — 29.35J calculated with the de- 
terminant method (det(VG'' — 1) = 0). In Fig. SI we plot 
the SPA approximation together with the exact results, 
showing that they are in good agreement everywhere ex- 
cept in the immediate neighborhood of the sharp reso- 
nance at [/ = f^ciR' ■where evidently the SPA fails. 

We now turn to the more interesting case of a quasi-lD 
two-body collision under transversal harmonic confine- 
ment, V{y) — ^y^- We solve Eq. (|3ip numerically and, 
doing so, we need an upper cut-off Nc in the number of 
single-particle trapped states. To obtain well-converged 
results, we calculate Xo.o = Uid/U using its Born series, 
Eq. ([5^ . to order 100 for different values for Nc, in the 
region where the Born series is well-defined. We choose 
Nc in such a way that the Born series offers converged 
results (not to be confused with the convergence of the 
Born series itself, which is granted) , and compare it with 
the exact numerical solution of Eq. (PT|) . The advan- 
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FIG. 5. Inverse tangent of effective interaction as function 
of the bare two-body interacting strength for a transversal 
harmonic potential with 0,1 J = 10"'' and K — Q. Red dots, 
blue squares and black triangles are calculated with different 
steps in \J j J to capture the width of the three resonances. 
Lines are guides for the eye. Insets: zoomed-in regions near 
the narrow CIRs, with axes labels the same as for main figure. 



tage of using the Born series for this purpose is that it 
is capable of handling much larger system sizes than the 
exact solution more efficiently when it converges. More- 
over, for very low values of ^2/ J, it is the only method - 
together with resummation techniques if necessary - that 
may overcome finite-size effects. In Fig. [5l we show the 
effective interaction C/id as a function of the bare interac- 
tion strength for a shallow harmonic trap (fi/ J — 10~^), 
which can be considered as almost separable into CM and 
relative coordinates [22|. The cut-off N^ — 41 is chosen 
since results are already well-converged. In this situation, 
we observe that a broad resonant profile is present, al- 
most unaltered by two very narrow CIRs that exist due to 
non-separability of the coordinates. Around these reso- 
nances, the effective interaction does cross zero and there- 
fore, in contrast to the separable case, true ID resonances 
are found. Since, of the three resonances observed, two 
are very narrow and do not seem to contribute much to 
the "overall" shape of the red curve in Fig. [5l the SPA 
approximation, Eq. (|33l) with U^f^ replaced by C/ciR' 
should give a good estimate of the broad resonance po- 
sition. Proceeding in the same way as for the two-state 



r(3) 



trap above, we find t^cm/"^ ~ —4.8 with an error of 
O(10~^), which indeed compares very favourably with 

the exact value of U^(^-^/J = -4.792. 

If we choose a higher trapping frequency, therefore 
making the separability assumption less rigorous, we ob- 
serve more narrow resonances, although some of them 
are much wider than in the almost separable case. In 
Fig. [HI we plot the effective interaction as U/J is varied 
for a trap strength of fl/J = 10"^ The cut-off N^ = 21 
already gives converged results. In this case the SPA 
approximation gives, for the broad resonance position a 

value of Uqi^/J w — 8.3 with an error ofO(10~^), while 
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FIG. 6. Inverse tangent of effective interaction as function of 
the bare two-body interacting strength for a transversal har- 
monic potential with fl/J = 10~^ and K = 0. Red dots, blue 
squares, black triangles and green diamonds are calculated 
with different steps in U/J to capture the width of the four 
resonances. Lines are guides for the eye. 



the exact position is U^^-^/J — —8.286. 

We note that a continuum of transversal states and fi- 
nite quasi-one-dimensional space for the two-body prob- 
lem are dealt with in the same way as in the potential 
scattering approach of Sec. |IT1 



IV. CONCLUSIONS 

In this paper, we have studied quasi-one-dimensional 
scattering in two dimensions with one of the directions 
tightly confined. We have used, instead of a continuum 
theory, a lattice model, and seen that single-particle po- 
tential scattering under transversal harmonic trapping is 
analogous to its continuum counterpart. We have then 
considered generalized situations in which the system is 
finite or the transversal confinement supports a contin- 
uous spectrum. For the two-body problem, we have de- 
rived the equation governing the scattering properties 
of the system, and this resembles a usual Lippmann- 
Schwinger equation which can be solved via standard 
methods from scattering theory. The center-of-mass and 
relative coordinates cannot be separated, and therefore 
our model accounts for the existence of more than one 
confinement-induced resonance. Moreover, zero effec- 
tive quasi-one-dimensional interaction strengths can be 
achieved, as opposed to the free-space, separable case. 
Our results are relevant for current on-going experimen- 
tal research on dilute ultracold gases in tight-binding, 
anisotropic optical lattices, where these systems are ex- 
tended over many lattice sites in the quasi-free direction. 

There are still many open problems in reduced- 
dimensional systems. Generalization and application of 
our results in a mean-field theory, e.g. in a transver- 
sal trap with two bound states, seems to be a natural 
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extension of the present work. We also note that an ex- 
act treatment of the bosonic and fermionic three-body 
problems may show modifications on the position and 
properties of two-body resonances, and may as well need 
genuine three-body effective interactions for their correct 
dimensional reduction 
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